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ON THE OSCILLATION OF THE LATERALLY 
HETEROGENEOUS EARTH, I 
Peter Musen 
Abstract 

The perturbative effects, as caused by lateral inhomogeneities in the Earth 
structure and by Coriolis force, contaminate the originally toroidal and spheroidal 
Earth’s oscillations, making them of mixed type. For this reason, in order to 
make the computation of the perturbations more uniform and homogeneous, we 
suggest the expansion of Earth's free oscillations into a series in terms of 
generalized harmonics familiar from the theory of angular momentum in quantum 
mechanics. Making use of Gibbsian symbolism and of some operators from 
the theory of angular momentum, we deduced the explicit expressions, in terms 
of generalized harmonics, of the perturbative terms in the differential equation 
of Earth’s free oscillations. We also obtained decomposition of the strain tensor 
in terms of canonical vectors. The integration problem for the cases of 
geophysical interest will be discussed in subsequent reports. 
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ON THE OSCILLATION OF THE LATERALLY HETEROGENEOUS EARTH, I 
INTRODUCTION 

In the present and forthcoming article we suggest an apparatus and scheme 
for the computation of the perturbed elastic oscillations of the laterally inhomo- 
geneous non-spherical Earth. The perturbative effects treated in this work are 
additive to the solution of the "main problem" as obtained by Takeuchi (1950) 
and by Alterman, Jarosch and Pekeris (1959), which represents the foundation 
of all modern theories of the long period oscillations of the elastic Earth. 

Methods being used to obtain the solution of the zero order are, more or less, 
a paraphrase of the Pekeris method. These authors have assumed that the 
density of the Earth’s material and elastic parameters depend only upon the 
distance of the particle from the Earth’s center and that the perturbative effect 
of the geostrophic force can be neglected. Under these suppositions the oscilla- 
tion can be conveniently decomposed into a sum of spheroidal and toroidal 
oscillations with the coefficients depending only upon the radius -vector. . 

The initial differential equation of the elastic oscillation can be split into 
disjoint final systems of ordinary differential equations for the coefficients. 

Each system must be integrated numerically using an adapted model of the Earth. 
It has now become evident, however, that the results of modern seismic obser- 
vations warrant the introduction of perturbative forces into the original Pekeris' 
theory. 
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The influence of lateral variations in density and elastic parameters (Toksoz 
and Anderson, 1966), (Dziewonski, 1970), (Lopatina and Ryaboy, 1971), and the 
perturbative effect of the geostrophic force, which causes the splitting of fre- 
quencies (Backus and Gilbert, 1961), (Pekeris, Alterman and Jarosch, 1961) 
are now being recognized in seismic records- The ellipticity correction, together 
with numerous coupling effects, became important and shall be considered. On 
occasion the Coriolis coupling between spheroidal and toroidal oscillations can 
exceed all other couplings (Luh, 1974). We must also consider the degeneracy 
and quasi-degeneracy of frequencies, familiar from quantum -mechanics. With 
these considerations the expansion of the perturbed elastic vibrations in terms 
of spheroidal and toroidal vector -harmonics can still be achieved (Arkani- 
Hamed, 1972), (Madariaga, 1972), (Luh, 1973), (Dahlen, 1968, 1969). 

It is doubtful, however, that the splitting of the displacement into toroidal 
and spheroidal components provides the best possible theoretical and compu- 
tational approach in the case of a rotating, laterally non-homogeneous non- 
spherical Earth when numerous couplings, degeneracies and quasi -degeneracies 
do occur. It seems that decomposition of oscillations into a sum of so-called 
"generalized spherical harmonics" can provide a better service, because it 
removes the discrimination between different kinds of oscillations and thus paves 
the way for more direct and homogeneous computational procedures- The 
generalized spherical harmonics are familiar from the theory of angular 
momentum in quantum- mechanics. There are numerous recursive relations 
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between them, and their products, which appear in the perturbing forces, can be 
conveniently expanded into Clebsch-Gordon series. 

It seems that Petrashen (1949) was the first to recognize the usefulness of 
the generalized spherical harmonics in the theory of elasticity. In recent 
years several works on the spectroscopy of the Earth were published which 
make use of the generalized spherical harmonics (Burridge, 1969), (Phinney 
and Burridge, 1973), (Smith, 1974) and the theory of group representations from 
quantum-mechanics. Formation of the scalar differential equations of Earth's 
perturbed oscillations is usually based on the application of covariant differentia- 
tion of vectors and tensors involved, using Einstein summation convention, some- 
times partially combined with Gibbsian symbolism (Backus, 1967). In the present 
article we suggest the use of Gibbsian vectorial and dyadic symbolism, and of 
operators familiar from the angular momentum theory in quantum mechanics 
to establish the differential equations of the perturbed elastic-gravitational 
vibrations of the Earth. The oscillations are decomposed along Petrashen's 
vectorial harmonics. 

The combined symbolism represents a fast and expedient geometrical way 
to establish the vectorial differential equations governing the oscillations and 
to convert them into the scalar equations. In the case of perturbed oscillations 
of the Earth these differential equations are no longer disjoint. The resulting 
scalar equations can be considered as a generalization of Pekeris equations and, 
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like them, can be integrated only numerically using an adapted Earth model. After 
the solution is obtained the oseilLations can again be easily split into spheroidal 
and toroidal parts. 

BASIC RELATIONS 

In this section we summarize some basic formulas to be used in the 
exposition. By applying the del-operator 


„ 3 1/3 3\ 

3r + r V* c)6 + sin 6 3tf>/ 

to the local unit vectors e 0 , , e f of the spherical system we deduce 

Ve ff = + ~ (- e^e r + cot e) , 

Vc ^> = " 7 (+ e # e r + e 4?d cot 6 )’ 


( 1 ) 


Ve r - + 7 (* e *e + ****>' 

and, as a consequence, by forming the scalars and vectors of the symbolic 
dyadics (2), we obtain: 


V*e„ = + 


cot 6 


V'e^O 


V • e 


r 



( 3 ) 
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„ 1 

V X = + - e^, 

Vx i (-e e + c r cot 6), 

V x e r = 0. 

Repeating the procedure, we have: 

V 2 e* = + 1 (-e*-2e r s in 6 cos $), 

r 2 sin 2 6 


and 




r 2 sin 2 6 


e <h> 


= v 

r* 


VV ■ (e 0 + e r sin 9 cos &) t 

r 2 sin 2 <9 

W ■ = 0 


(4) 


(5) 


( 6 ) 


W« e r = -■ 


We obtain a more compact system of formulas if we decompose vectors and 
tensors which appear in the present exposition not along eg, , e f , but along 
the isotropic canonical vectors 
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and along 


*- = + — (+ e s “ ie s) 
/2 


e 4- = (+e fl + ie^) 


+ “ /2 


e 


o 


e 

r 


( 7 ) 

(8) 

( 9 ) 


It is convenient in the frame work of the present exposition to define the 
the products of (7) - (9), scalar, vectorial and dyadic, in accordance with the 
standard rules of classical vector algebra, considering the complex coordinates 
as scalars. 

From the orthonormality of e 9 , e^, e r we deduce: 
e_ • e_ = e + • e + = e_ • e 0 = e + • e 0 = 0, 

, _ i (10) 

Cq*6q — + 1, 

e + x e_ = + ie 0 , e Q x e_ - +■ ie_, e + x e Q = + ie + . 

Elimination of e 0 , e^, e f from the idemf actor 

I = e 0 e g + + e r e f 

in favor of canonical vectors leads to the representation 

I -- e 0 e 0 — e + e_ - e„e + 

which is useful in performing the rotational transformation of vectors and tensors. 
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From (2) - (6) we deduce a set of basic relations: 


Ve_ = + - 
r 


- e_e 0 — L- (e_ + e + ) e_ cot 6 

V 2 


Ve + x + -L 
r 


- e +«o +-4r ( e - + * + ) e + cot # 

/2 


(ID 


v e 0 = -- (e_e + + e + e_), 


and, by forming the scalars and the vectors of the symbolic dyadics (11) and 
taking (10) into account, 


V • e_ = + 


cot 6 
r v’2 


V-e + x- 


cot 6 

T/F 


( 12 ) 


V-e n = 0, 


^ x e. = + - ( K. cot 6 

r \ /2 


Vxe + = + — (-e + cot 9 

/2 


(13) 


V x e 0 = 0, 
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V 2 e_ (e_ + e Q y/2 sin 0 cos 8), 

r 2 sin 2 8 


V 2 e + = (e+ - e 0 V2 sin 8 cos 6 ), (14) 

r 2 sin 2 8 



W 


* e_ = - 


r 2 s in 2 8 L 


1 e o 

+ (e - e + ) + - sin 8 cos 9 

2 ~ ' v2 


VV - = + 


1 


r 2 sin 2 8 L 


1 e o 

+ -(e -e + ) + ~sin(9cos0 

2 /2 


(15) 


VV- 


2 


e 


o* 


The set 

Yt - e 
~ + •». 

Y ^m = e - Y «>. +1 (*.*>’ 

v d _ ~ y ^ 

, '£ m " e 0 * m, O(0,^>’ 


(16) 


s ? j) 

where Y , Y Y „ are particular elements of the set of generalized 

m , 1 m , + 1 m , 0 

'i 

spherical harmonics Y^ , familiar from angular momentum theory, and they 
constitute the basis of the expansion of the elastic displacement of the Earth 
(Phinney and Burridge, 1973), (Petrashen, 1949). We have 
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(cos 8), 

m , n mn v ' * 

where P mn (cos 0) are the generalized associate Legendre functions (Edmonds, 
1960), (Vilenkin, 1965). 


The selection of (16) as a basis for the expansion removes the discrimination 
between spheroidal and toroidal oscillations and makes the computation of 
coupling effects more uniform and homogeneous. After the solution is com- 
pleted we can again decompose it into spheroidal and toroidal modes. When no 
ambiguity results we will omit, for the sake of brevity, the indices and m 
in (16). 

The reduction of the vectorial differential equation governing the oscillations 
to the scalar equations can be simplified by making use of the operators 


H + = e*^ 


+ cot 8 


JL 

Bl/i 


1 3 . _a\ 

sin 6 B<£ + 1 B 6)' 


H_ = e +i ^ (- cot 8 — 

\ cM 


1 


cty sin 8 


3 . B> 

— + i — I 

B <p ~d8) 


(17) 


of infinitesimal rotations and by introducing the generalized del-operator: 


V = e 0 — + (e"^e + H_ - e +i ^ e.H*), (18) 

3r r/2 

where i p, 4> and 8 are Euler angles. 


By applying the operator (18) to the elements 

m,n jn,n 

of the irreducible representation of the rotation group and taking into account 
(Edmonds, 1960), (Vilenkin, 1965) 
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where 


and 


0 £ n = [(■£ + n) ('t - n + 1)3 1/2 , 


Ve -in ^ = - _£L- (e + + e.) cot e e ~ in 'f', 
rv"2 

we deduce after easy transformations: 


V = + — - — (etc e Y^ +1 - Q-o e + Y_ .) 

m,n T ^ 'V ,n+ l ~ 01 > n + * / t/,n + m « n * 

lV2 


+ + e + ) Y cot 6. 

r/2 


(19) 


Taking (10) into account, we have 


* x VY * — — — — (cl o c Y . . + 

'0 * v 1 m, n v l-.n+l m,n+l 




ao e + Y )+-iiL(e 
H-sn + 01,11 1 r /2 


«+> Y ln COt 8 

( 20 ) 


In particular 


ia. 


VY r + 

V 1 m, 0 + 


«,0 


:v£ 


<*i. - Y l->’ 


and 


( 21 ) 


^ n.fl ^y- Y + 

e " X VY - - <-*l m + Y 4n->- 


0 * v ■ l m,0 


r ^2 


(22) 


The left-hand sides of (21)-(22) are the spheroidal and toroidal harmonics, 
respectively. Thus, the generalized spherical harmonics Y~, Y + represent 
the linear combinations of the standard spheroidal and toroidal harmonics. 


10 
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Making use of (19) and (11) we deduce: 


VYt = + ~ (a, e_Y^ „ - a» e.Y^ ,) e. 

r/J ,0 "t.-l + m > 2/ + 


i 

" 7 Y m, -l e + e 0’ 


= + ^f (a l2 C ' Y M2 ~ a 4 ,!*+*«.<>> «- 


1 

“ T m, +1 e - e 0’ 


(23) 


VYS = - lY^ n (e_e + + e + e.) 

'v m j-m,0 + + 7 


+ -75 


r/2 


and, as a consequence, taking (10) into account: 


Vx Yi = +-L2 y $ -Ay), 

1/111 rv^ ^ r 


Cl j? 

V7 V ** ^ i 0 *| 7 0 1 17 

v X )[ — + '■ 1 IP + — \ o » 

7T 4s m x* ^ m 


(24) 


r/2 


vx y I = + !£lL(y ; +y; ), 

'I'm r/2" 
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and 


v-y; =-v*Yp 

/ T/ m 'v m 


lCt ^ , 0 y't V . Y 0 _ + 2_ Y '£ 

~ 'W 


ri/T 


i,0 


r m , 0 ’ 


(25) 


DECOMPOSITION OF THE STRAIN TENSOR 

In this section we deduce the canonical dyadic decomposition of the strain 


tensor 


e " 2 ( uV + Vu ) 


(26) 


assuming the canonical decomposition 


u = U° (r) Y° + U + (r) Y + + lT(r) Y 


(27) 


of the displacement. 
We have 


Vu ~ e 0 

0 \ dr dr dr 


+ (U°VY° + U*VY + + U'VY") 


(28) 


Taking (23) into account and rearranging the terms in two different manners, we 
obtain the following two canonical dyadic representations of the gradient of the 
displacement: 


Vu = 


HI1 U 10, n 

-^-eoYo + -pU°(e.Y +1 -e + Y_ 1 ) 

rv'2 


U e Y - U “ e Y 
e + i_j - — - e_Y +1 


(29) 
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Eh! e Y - H! e Y 

tdT 0 “ J ~ 0 


+ -i— U + (a 0 e.Y 0 - a - 1 e +Y_ 2 ) e+ 
r/2 


dU' 


L dr 


U u 

e oY +1 - 7- « + Y 0 


and 


+ _U-(a 2 e.Y +2 
r/2 


- a i e +Y 0 )j 


r, _ /dU° v dU* v dU" v 

+ C ° Vd7 e ° Y ° + dr C+Y_1 + d7" e_Y+1 , 


+ e.i 


; u “ ( + ? e - Y ° + ^ •» Y ->) 


( 30 ) 


- U + ( + — e + Y_, + - e 0 Y_ 

V r/2 


la o _ 
~U e.Y 0 

r/2 


+ e_ 


/ i ia n 


ia A 

+ U + -ie + Y 0 

r/2 


+ U (+ e.Y +2 -ie 0 Y +1 

V r/2 r t 
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Combining (29) with (30) transposed, we obtain the canonical decomposition 
of the strain tensor: 

e = ( e °° e o Y o + e+ ° e + Y -i + e "° e - Y +1 ) e 0 

(31) 

+ + e ++ e + Y _ 2 + e _+ e_Y 0 ) e + 

+ ( e ° e 0 Y +l + € P+ Y 0 + € e - Y +2^ ®-’ 

where 

,oo _ dU° 

' d7’ 


ia. 

e ++ = 1 U + , 




la 2 - 
e = + — ±U , 




e +o _ e o+ _ I dU + U* la o yo 


2 \dr r r /2 


(32) 


e“° = e°~ = I ( SI - HI + u° 


2 \ di 


r/2 


T T o ia n 

e _+ = e + ~ = - — + — (U + - U‘) 

r 2 r/2 

By forming the vector of the dyadic (29), i.e. replacing all dyadic products 
by the vectorial ones, and taking (10) into account, we deduce 


Vx u = Z°Y° + Z + Y + + Z~Y“, 


14 



17 


where 


z° = +~L(U + + U“) 


z + = + -iu°.i 


dU + U H 


Z~ = + — — —U° + i dU + 


dU~ IT 
dr + r 


Similarly, by forming the scalar of (29), i.e. replacing the dyadic products by the 
scalar ones, we have 


where 


V-u = X°Y 0 , 


X«(u) = ^! + ^! + ^(U--U + ). 
dr r r/2 


Repeating the process for V x u , we obtain: 


where 


V x V x u = W°Y 0 + W + Y + + IT, 


w° = + 77f (z+ + z ")> 


W + - + Z c - i m + 

r/2 ' dr r 
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From (34) and (21) we deduce 

W-u = — Y° + X°0T -Y*), (38) 

d ' r/2 

We shall need (29) - (38) in the transformation of the differential equation of 
oscillations of the Earth. In this process we shall also make use of the Clebsch- 
Gordon expansion 

K*pI 

Y l Y qs = Xl C( ' t ’ P ’ ^ m ’ q) C ^’ P ’ V: n ’ S) Y »+<J.“ + - (39) 

V = \-t ~ p | 

where the symbols C(-£, p, v ; a,,/3) designate Clebsch-Gordan coefficients. They 
can be computed directly on an electronic computer or can be obtained from 
existing tables. 

Every scalar in the present theory can be represented as a series in har- 
monies Y m 0 and every vector as a series in Y^ , Y£ , Y^ (-1 = 1, 2 , 3 , . . . ; 
m = The coefficients in the expansion depend only on r. 

DIFFERENTIAL EQUATION OF THE ELASTIC OSCILLATIONS OF THE EARTH 
If we assume hydrostatic equilibrium, but permit lateral inhomogeneities in 
density and in elastic parameters, then the differential equation of the elastic 
oscillations of the Earth takes the form: 
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p — + 2pfikx = V'OT- v-(pu) vv 

Bt 2 

+ V(pu * VV) + pVi/<, (40) 

where p is the undisturbed density, -V - the interior gravitational potential (it 
includes the potential of the centrifugal force). Cl - the angular velocity of 
rotation of the Earth, k - the unit vector alongthe polar axis, - ^ - the increment 
in internal geopotential due to redistribution of masses, and a -the stress tensor, 

cr = nV'n + p(uV + Vu). (41) 

We assume that the density p and the elastic parameters X, p. can be represented 
in the form: 

p = p 0 ( r ) + S P> 

(42) 

\(r) + S\, 

M = ^o( r > + S * a > 

where the perturbative terms 2p represent the lateral deviations from 

the spherically symmetric mean values p 0 , ^ 0 , p 0 . We assume the existence 
of the expansions 

+ CO q 35 + p 

ip= L L p ^ r)y ^- 

P’l q = ~ P 

+ co q—+p 
p = 1 q*-p 
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+ oo q — + p 
P =1 q=~P 

(p= 1, 2,...; q = - p, ..., + p) 

in terms of associated Legendre functions For the geopotential inside a 
given domain of integration we assume an expansion of the form 

V - V 0 (r) + SV, 


5V = 21 'WOY*, 0 - <43 

& 

With the domains of convergence determined by the shape and position of the 
surfaces of discontinuity inside the Earth. In addition to (42) we found it useful 
to introduce the combination: 


fi = /3 0 + S/S = X + 2fx = /3 0 (r) + 


+ co q “ + p 



pq 


< r > YP q , 


(44) 


A) - \ + 2 ' t V ^pq = \>q + 2 ^pq- 

Substituting (41) into (40) and taking the Fourier transform (or searching for a 
periodic solution) we obtain 


pi-a?! + 2ifL>k x I) ■ u = /3W- u - /uV x V x u 

+ 2e ■ Vjt + (V- u ) V\ - V- (pu) W + V(you ■ VV) + p V</r, 


(45) 


where u now designates the transform of the displacement. 
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We assume that the higher order terms 

SpVSV and (VS p) (VSV) 

are negligible. Under this assumption the equation (45) takes the form: - 

p (- &j 2 1 + 2i<nfik x I)'u = /3VV*u- p.V x V x u 

+ 2e • Vp + (V* u ) V\ - V- (p 0 u) VV 0 + V(p 0 u • VV 0 ) + pV0 ' 

- [V* (uSp) VV 0 + V- (p 0 u) VSV - V(Spu • VV 0 ) - V(p„u * VSV)] . ( 45 ’> 

The last four terms are perturbative. They represent the effects of the lateral 
inhomogeneities in density and of the deviation of the figure of the Earth from a 
sphere. 


UNPERTURBED PROBLEM 

The solution of zeroth order is obtained under the assumptions that the 
density, elastic parameters and the interior geopotential are spherically sym- 
metric functions, that the perturbative effect of the geostrophic force is negligible 
and the oscillations are purely elastic. For the solution of the zeroth order we 
assume an expansion of the form: 


u = 



(4&) 


where 


U 'tm = U L (r ’ Y <L + U lm (r ’ Y L + U ^ m (r> 


(47) 
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and the eigenfrequency o> is determined from the boundary conditions on the 
surfaces of discontinuity, assumed to be welded together (if the materials are 
solid). The radial functions u£ , U|> + , satisfy a set of disjoint ordinary 
differential equations. For each u„ they can be integrated separately. For the 

“T/ m 

zeroth order solution of the form (47) we have for the particular terms in the 
differential equation (45), taking (31) - (38) into consideration (we repeat here 
the complete set of formulas given previously) , 


where 


/3W 'u = /3 0 


jvO i (X n 

y° + — X °(Y~ 

_ dr r/2 


-Y + )| 


(48) 


dU° 

d7 


x° = ^ + - u° +■ 

r 


ia. 


r/2 


- (IT - U + ), 


(49) 


/xV x V x u = /Xq 


+ — <Z + + Z") Y° 
. r/5 



(50) 


+ 






where 


Z° = + — (U + +IT), 
r/2 


Z + 

z~ 



( 51 ) 
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where 


2e • V/x = 2 _JJ (e 00 ¥° + e 0+ Y + + e°~Y~'), 

dr 


£ oo = dU_° 


.0+ = i + _^V 


2 Vdr r ri /2 


£ o- 1 te.tr^uoy 

2 \ dr r iV2 / 


(V-u) V\= — 9 X°Y°, 
dr 


V.(p„)VV = ^(^ U° + P 0 X°J Y°, 


and, taking (21) into consideration, we have: 


V(p 0 „.vv) = £ Lu« ^)y° -£ku* £) cv- -»*)• 


For the increment of the interior geopotential we have 


= P(r) Y 0 , 


assuming the Earth to be a selfgravitating system. From (57) we have: 


dP v0 ia o 


V</; = ^ Y° +— P(Y~ - Y + ), 
dr r/2 
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and, taking 


a n 

V 2 Y - Y 

v x 0 i 0 

r 1 


and (25) into account, the Poisson equation 


V 2 ip = + 47 tGV • (jOu ) 


(59) 


becomes 


£E + 2^-!£p = 4,g(u»^^ 0 x»), 

dr 2 r dr r 2 \ dr / 


(60) 


the same as in the Pekeris work. For the sake of brevity we systematically 
omitted the indices -£ and m in the equations (48) - (60). 

Substituting (48), (50) or (36), (52), (55), (56) and (58) into (45 r ) and keeping 
only the terms of the zeroth order we obtain the differential equation: 


p^cJu = 




dfx 0 


dK 


dr 


dr 


dr 


f 2 dv o d2v o 


P 0 \~ d 1 U ° ~ — F “^o( u ‘’* u+ > + /°o^7 

°Vr dr dr 2 / rvf dr dr 


ia A 

+ -4 ( / 3 0 X° +Po P) (Y--Y + ) 
r/2 . 

+ 2 ^ e°* - M 0 W" Y' + (2 ^-° e 0+ - M 0 W + ) Y + , ( 61 ) 
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or , in scalar form 




dr 


(62) 


/ 2 <*V 0 d 2 V 0 \ ia 0 dV 0 dP 

'H-ff-W 0 -^^ (u - u > + p »d7' 


-p 0 a?\J~ = + K° -yu. 0 r + 2 -£? e 0_ , 

- / o 0 « 2 U + = - K° - m 0 W + + 2 ~ e°\ 

dr 


(63) 


where 

Kl> = + ^ f° X ° + p ° u ° 37 + PoP ) <64) 

and X° is given by (35), W°, W", W + by (37) and (33), and e 00 , e°~ , e 0+ by 
(32). 

PERTURBATIVE TERMS 

In this section we give the typical perturbative terms in the differential 

equation of the elastic oscillations (45'). These typical terms carry the influence 

on u ; of the coupling effects between p , K , u and u , and of the geo- 
'im p q pq 'pq sk 

strophic force. In developing the formulas we assume that the square of the 
displacement (and frequency) can be obtained by means of the Schrodinger- 
Rayleigh technique as a linear combination of Petrashen’s harmonics. 
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From (38) and (44): 


F, = 0VV- u 


+ i; 




(l ; m k q’ °) = + C< ' S ’ P ’ k ’ C ^ S ’ P ’ Sk +£ i- r 


pq dr 


m k q ’ 


itt a 

= + C(s, p, -£; k, q) C(s, p, + 1 , 0 ) S k+qim / 3 pc( — — X^ k , 

rv 2 


(lj mkq’ + ) = ~ k, c i) C(s, p, - 1 , 0 ) S k+1 m /? p 


so v o 


and from (36) and (43): 


F, = /LiV x V x u 


{ s p, 


- + 2; 


m k q 


• °)n.* ( 2; 


Z s p. 

m k q 1 


*<.♦ kM i( ) n.- 


where 


( 2; t k q ; °) = + C(S ’ P ’ l '' k ’ q) C(S ’ P ’ * 5 °’ 0) ■W-''*"** 

( 2; m k q’’ -) = + C(S> P ’ l ' k ’ q) C(S ’ P ’ l ' + 1 ’ 0) WnAq^k’ 
( 2; m k q 1 + ) = + C ( S> P ’ ^ k ' Q ) C ^ S ’ P ’ ^ _1 ’ ^ S k+q,m^pq W sk ' 


From (21), (31) and (43): 


F 3 = 26-V^ 


= 3; 


m k q 


I s p. 
m k q ’ 


*<.♦(* i w 
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where 


3; ^ f p ; 0 

m k q 


+ 2C(s, p, ■£; k, q) S 


k +q t m 


+ C(s, p, h 0, 0) 


^p q [ 6 “k C ( S - P’ *'• "I- + 1) - € sk" C ( S > P- +*> - 1 )] 


. . -£ s p 

’ m k q 


; = + 2C(s, p, Z] k, q) S 


k+q , m 


+ £ rtC(s, P, V, + l, 0) 


A U. q 

“ ^pq ^ sk ^ ’ P* + ^)~^ s kC(s, P» ^ + 2 , 

r/2 


( 3; fkq 1+ ) = + 2C(S ’ P - * k ' q)S k*0.. 


pq - +o 


e sfeC( s > p, - i, 0) 


“-j^Aipq [< k + C(s, p, -2, + l)-e£C(s, p, •€; 0, _l)]i, 


where, in accordance with (32): 


UL VJ , 

: 00 _ sk 


tr + 


€ sk = + 
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. +o _ ,-0+ _ 

sk sk 


€ -o _ € o- . 

sk ~ sk “ + 


1 j 

f d «:k 

_^k_ 

ia sO 

2 

^ dr 

r 

r/2 

1 1 

( du ; t 

u; k 

Ko 

2 ' 

I'dr 

. . 

r 

r/2 


u 


sk / * 


u 


sk 


+ - e + " 

sk sk 


»!k 


- + — ~ (u* k - u ; k ). 

r 2rv'2 


From (21), (34) and (43): 


F 4 = (V-u) VX 

= (^k£°)n. + (r i t s * +) n.- 


where 


dX 

o . p°i 


4; m k q : °) ~ + C ( S ’ P ’ k ’ q ' ) C< - S> p! ^ °’ °' ) S k+ q ,m X S k "57“ ’ 

4; Ik a : “) = + C ( S) p ’ k ’ q ' ) C(s ’ p ’ ^ °’ +1) Sk+ 'i- m “ 

v m k q / rv ^2 


0 

sk 


( 4: m k q ; + ) •- <=<*• »• * k. <0 c&, P, t: 0 . - 1 ) s l4 ,.. ^ k„,x» 

and again from (21), (34) and (43): 

F,=v-(„sp)vv 0 = (s;^q ; 0 ) ¥ L' 


where 


ST i ’ p ; 0 
m k q 


= C(s, p, l\ k, q) S k+q>n 


dv 0 

dr 


'dp, 


+ P 00 X° C(s, p, 0, 0) 


- 57 - ~sk T ' pq sk 


ia 


4 _^p pq [u; k C(s, p s t; +1, -l)-u; k c<s, P, 1-, -1, +1)] 
rv ^2 




(71) 


(72) 


(73) 


(74) 


(75) 


26 



29 


F 6 = Vv (P 0 U ) VSV 


= 6; 


s p 
mkq 


; °) ( 6l «Jq ; -) Y i. + ( 6 -fkq ; 


; +J y; , 

'T/m 


where 


( 6; mk? ; °) = + C ( s > P- l ' k ’ q > C ( s ’ P’ ^ °- °) S k+q ,< 


dp, 


• 


dV 


6 ; 


. 'C s p. 


mkq’ 


= + C(s, p, 4; k, q) C(s, p, l\ 0, + 1) 5 k+q!n U° k + P 0 X S °J < 76 > 


ia 


P o 


V 

- — pq' 

r/2 


( 6 ’ m k q’ + ) C ^ s ’ p7 k ’ P? S k+l (“37 + ^° X sk) 


ia 


p0 


v ~ 


F 7 = V(5/3u-VV 0 ) 

= (7; 1 f P ; 0) Y? + (7 ; ^ f p ; -) Y» + (7; * ? p ; +) 

\ mkq / -vm \ mkq / -tm \ mkq / 


Y L> 


where 


( 7; m k q ; °) = + C ( s ’ p ’ k ’ q ) C ( s ’ p > ^ °> °) s 


k+q , m 


If u» S 

dr V°P<f * k dr/ 


( 7 ’ m k q* ~) = " ( 7: mkq’ + ) = + C(s ’ p ’ ^ k ’ q) C < s ’ p ’ °’ °> S k+q,m ^ 


ia 


r/2 


f, u° dV »' 

«- VPQ S 


M 3k dr; 


27 



30 


where 


and 


8; 


F 8 =V(p 0 u-V8V) 

= ( 8; mkv °)n.+ (* 


Z S P. 

m k q ’ 


n,* i e ^kq^ *) »*,• 


ft. I s p . n \ _ d fo. ^ s pi 
, 8, m k q’ / dr ’ no k qj ’ 


m k q 


1 s p. ,'i - , ‘^o r 8; *- * pi 


; -) =- jg; "“P; J =+ _i» Ifc, . 

/ V ra k q / rv /j |_ rnkq 


> 

ij 


"£ s p 
rnkq 


f dV 

= p 0 C(s, P> -t; k, q) S k+q m iC(s, p, 4; 0, 0) U° k + 


(78) 


(79) 


la 


+ — ^ ^o v P J c ( s . Pi h +li ~l)U; k - C(s, p, Z; -1, + 1)U* ] 


r/2 


In order to obtain the representation of a typical term in the expansion of the 
geostrophic force 


iF 6 = 2 p 0 io;nk x u 

We substitute into the last expression 


k = - e d sin 6 + e f cos 9 


or 



(e + - e_) sin 6 + e 0 cos 9 


and 


u - u 


4-i. a + a Z, m + tt Z+i, t 
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Taking 


pi _ pi pi _ pi _ 

o,+i r + i,o r o,-i r -i,o 


i s in 9 

v2 


Pq 0 = + cos 6 


(80) 


into consideration and making use of the Clebsch-Gordan expansion, we obtain 
for the typical term in F 9 : 

f - = (* 1 °) n. * ( 9; i- -) - (* m : + ) n.- 

where 

(9; l m , 0) = C(4-l, l,4;m, 0) [C(4-l, 1,4; -1, +1) U|. 1 m +C(4-l, 1 ,4; +1, -l)U^_ x J 
+ C(4, 1, -t; m. 0) [C(4, 1, -t; — 1, + 1)U| ra + C(4, 1,4; +1,-1) iTj J 
+ C(4+1, 1, 4;m, 0) [C(4+l, 1,4;-1, +1)UJ +1 m + C(4 + 1, 1, 4; +1, -1) U£ n J, 

( 9; m : ") =C (' e - 1 ’ °> EC(4-1, 1,4; 0, +1)U|_ 1 m + iC(4-l, 1,4;+1, 0)1^ J 

+ C(4, 1, 4; m, 0) [C(4, 1, 4; 0, +1)U| ^ + iC(4, 1, 4; +1, 0) J (81) 

+ C('£ + 1, 1) m, 0) [C(^ + 1, 1> 'fc r 0, + 1) + iC( r t + 1, 1 , 't; + 1, 0) Uj> ] 

( 9; m’ + ) = 1 - ^ m ' °) tC(4- 1, 1, 4; 0, -l)U|_ i m - iC(4- 1, 1, 4; -1, 0)U|_ j J 

+ C(4, 1, 4; m, 0) [C(4, 1, 4; 0, -1)1^ - iC(4, 1, 4; -1, 0) u| J 

+ C(4+ 1, 1, 4; m, 0) [C(4 + l, 1, 4; 0, -l)u£ m - iC(4 + l, 1, 4; -1, 0)U| +1 J 


Finally, for the perturbative term 


F , 0 = w 2 Spu 
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we obtain: 

f 10 = (io; 

where 


lsp t 
mkq 


; °) n. 




K:;=*) 


v + 


mkq’ a, oi 0 pq^'( s ’ P> k> c l)C(s, p, 'C; 0, 0) S k+q m , 

( 10; mkq ; ~) =a ^ P <, C < s > P- k k. q)C(s, p, 4; +1, 0) 8 k+qn)1 

( 10; m k q ’ +) = ^ Pq C ( s > P- k > q) c ( s ’ P. k - 1 - °) 8 k+ „, m 


(82) 


In the group of formulas (66), (68), (70), (73), (75), (77), (79), (83) the selection 
rules are 

js -p| f> ^ S + p 


k + q = m. 


and 


|k| S p, |m| 6 l. 


CONCLUSION 

The proper understanding of the mechanics of free oscillations of the Earth 
permits us to improve our knowledge about the internal structure of the Earth. 
Free oscillations are normally expanded into a series of vectorial and toroidal 
harmonics, assuming spherical symmetry of the Earth. However, the perturbative 
effects of lateral inhoraogeneities, of the Coriolis force, and of couplings between 
the oscillations are now also being considered in seismology. They cause the 
splitting of frequencies and "contaminate’' the originally pure toroidal or pure 
spherical oscillations. The originally toroidal oscillations aquire perturbative 
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spheroidal components and vice versa. For example, the Coriolis force adds 
radial components to the toroidal ones (Mac Donald and Ness, 1961). Thus we 
can talk only about oscillations predominantly spheroidal or predominantly 
toroidal. For this reason, and to make to computations of the perturbative effects 
more uniform and homogeneous, we suggest in the present work the expansion 
of perturbative effects in the free oscillations of the Earth in terms of vectorial 
generalized harmonics familiar from quantum mechanics. We give the explicite 
form of those perturbative terms in the differential equation of the Earth’s 
oscillations which are caused by the lateral inhomogeneities and by the Coriolis 
force. The problem of integration for the cases of geophysical interest and of 
degeneracy we shall treat in subsequent reports. 
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